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ABSTRACT 


A quasi-geostrophic omega equation for large-scale vertical veloci- 
ties in the atmosphere is solved at the earth's surface, 850 mb, 700 mb, 
500 mb, and 300 mb for a grid of 1977 points per level covering the major 
portion of the Northern Hemisphere. The effects of terrain and surface 
frictional stress on the vertical velocity at the lower boundary are 
investigated in some detail and applied at the earth's surface. The 
three-dimensional integration of the resulting equations excludes areas 
"inside" the terrain. The mathematical model is programmed for the 
Control Data Corporation 1604 digital computer to utilize the operational 
output of machine-processed data and analyses produced by the Fleet 
Numerical Weather Facility (FNWF), Monterey, California. The program 
allows a wide variety of modes of operation and yields a solution for W 
in five minutes. 

An unusual weather situation in January, 1962 is investigated and 
the results, though tentative, show that the correlation between computed 
vertical velocities and the observed weather is excellent. 

The writers wish to express their appreciation to Professor George J. 
Haltiner of the U. S. Naval Postgraduate School for his guidance, contri- 
butions and assistance in this investigation. 

Appreciation is expressed for the aid and cooperation of personnel 
at FNWF, especially to Mr. Milton H. Reese, who designed and programmed 
the initial phases of the machine computations. Gratitude is expressed 


to Mrs. Lillian F. Paine for her patience in typing this thesis. 
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LIST OF SYMBOLS 


The vertical velocity of air in pressure coordinates 

A subscript that specifies the value of a parameter at the 
lower boundary 

The frictionally induced W at the top of the friction layer 
The potential temperature 


The vertical component of the mesh Kg ); K increases 


with increasing pressure 

The universal gas constant 

The coriolis parameter, 22 sing , where (2X is angular 
speed of earth and P is the latitude 

The value of the coriolis parameter at the pole 

The mean value of the coriolis parameter 

The upward component of the apparent gravitational acceleration 
The relative vorticity 

The absolute vorticity 

The horizontal Laplacian operator on a constant-pressure 
surface 


The del operator on a constant-pressure surface 


The horizontal Jacobian operator 


) = 9A 0B _ 3B OA 
J (A,B) dx dy Ox Oy 


The height in centimeters 
A subscript that specifies the value of a parameter at the 
terrain 


The specific heat of dry air 





The specific temperature anomaly 
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Used as a Subscript to specify difference of a parameter 
between two heights 

The pressure altitude in centimeters 

The mesh dimension, 384 km 

The altimeter correction, DS LES 

The surface frictional stress 

The map factor for a polar stereographic map projection true 
at 60N 

A subscript that specifies the value of a parameter at a lower 
(L)/upper (U) height 

A caret specifies that the quantity has been scaled into a 


fractional quantity 
Te 
Q dp 


The atmospheric pressure in millibars 


The static stability, JO = - 


Finite difference operators 





is Introduction 

The value of vertical velocities in the atmosphere for understanding 
the basic weather phenomena associated with the large-scale features is 
gradually becoming more fully appreciated by synoptic forecasters [1]. 

The purpose of this investigation is to provide a means of computing these 
vertical velocity fields from routine weather charts and to examine the 
Structure of these fields. Among the various methods for such determina- 
tions, the omega equation, derivable from the vorticity equation and the 
first law of thermodynamics, appears to afford the best means of computing 
large-scale vertical velocities from synoptic meteorological observations. 
Thus, the computer program developed for this omega equation is adaptable 
to operational uSe in a numerical forecasting center as well as for basic 
research. 

The data used in this program are produced by the Fleet Numerical 
Weather Facility (FNWF), Monterey, California, essentially as follows [2]: 

(a) Northern Hemisphere weather data are processed automatically by 
digital computer to yield reliable listings of surface and upper-air data. 

(b) A Northern Hemisphere surface map is analyzed by digital computer 
and 1000-mb height and temperature charts are produced. 

(c) The 850-mb, 700-mb, 500-mb and 300-mb height and temperature 
charts are then constructed with hydrostatic consistency by building 
upward from the 1000-mb level. Winds and off-time reports are also 
included. 

(d) These surface and upper-air charts cover the entire Northern 
Hemisphere to the equator, and are produced at 0000 GCT and 1200 GCT 
daily. Further details are available at the Fleet Numerical Weather 


Facility. 





ae Method of Solution 

On the basis of the available data described above, it was decided 
to arrange the computational form of the omega equation so that only the 
data at a given constant-pressure level would be needed to calculate the 
parameters required in the iterative Scheme for that level. This necessi- 
tated the use of unequally spaced pressure intervals in computing 
vertical derivatives, but gave the model some flexibility if additional 
levels became available or were changed. 

The omega equation is derived from the vorticity equation (1), and 
the thermal equation (2) for adiabatic conditions by eliminating the 


local time derivatives jointly between the two equations. 
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where the various symbols are defined in Table of Symbols, page vy. 
The geostrophic approximation for vorticity has been made in the first 


term of equation (1). Operating on (1) with ee andion (2). With 


(V7 } and adding results in a form of the omega equation 
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The terms involving the vertical advection of vorticity and the 
twisting terms are usually of opposite sign and, moreover, small in com 
parison with the others remaining, and will be neglected as a first 


approximation. Of the terms which arise from the horizontal derivatives 


Ot GC, 
9 / 2. zs : 

MICO) =TVo+woVG+eva-vo (4) 
the last two on the right are also considered to be small in comparison 
with the first and will be neglected. The neglected terms may be 
incorporated along the lines used in Stuart [3], to afford a solution 
to the more complete equation (3). It is to be noted that 7 will be 
allowed to vary in the horizontal as well as with pressure in this scheme. 
To evaluate partially the effect of the variable &@, the computer program 
has been designed to utilize alternately horizontal and vertical means of 
the static stability as well. 

Neglecting the last two terms in (4) and substituting in (3) gives 


the reduced omega equation (5) 


2 2 
2. tT Bays Zo ee. ; Ve he 
o Vint Bln 2% _ PYyevT)- PL 2 (Vevz) =O (5) 
R op RK op 
2 

where, in the second term, fr has been replaced by the to meet certain 
consistency requirements as described in [4] and [5]. Again, to evaluate 

ca 
the effect of this term, the computer program has both versions Ty and T 
incorporated. 


The third and fourth terms of (5) can be expanded using the geo- 


strophic vorticity = g wz : 


i 





the hydrostatic equation, and the equation of state: 


Iv (V- -- ie a v7)| = 


9 2 ae a 2 
= |-v 5(2,7)-E5(50) +3(Z,V 7] 
Substituting this expression into (5) and rearranging gives the computa- 


tional form of the reduced omega equation: 


Vw+ Pla = = Slvr Zt) = (? \VT)- 17,0) (6) 
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The basic parameters involved above are Zz and ale - At a given 
level, the right side of (6) can be calculated from data at that level 
only. The last term on the right in (6) is generally considered to be 
of major importance and is used in the approximation made by Wiin-Nielsen 
[6]. 

The grid used is the Fleet Numerical Weather Facility Northern Hemi- 
sphere octagon of 1977 points centered on the pole and extending to about 
latitude 10N. 

The omega equation (6) is to be solved at the four levels of 850 mb, 
700 mb, 500 mb and 300 mb, with the lower boundary being the surface of 
the earth. The upper boundary is 100 mb where it is assumed that W=O , 
and the lateral boundaries are near 10N where it is also assumed that 
n=O . 

To solve (6) numerically, the following finite-difference approxima- 


tions, based on the adjoining five-point array, were uSed. 
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The portion of the term involving P Vom 1S approximated [4] by a finite- 
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difference analog for unequal intervals 
: S 
Peet. =P | 2 Aw (8) 
Z 
OP Ayp Ap 


where Ae = W K-| ee AxP LO + AvP Wee 
Dip Ay, p 


and AP » ALP and Ap are defined in figure l. 
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The static stability 


8 ed 9 (9) 


is equal to the difference between the dry adiabatic lapse rate 
bp = RT/p Cp in pressure coordinates and the observed lapse 
rate Yp = OT/ 2 P . This form of the static stability parameter 
corresponds to that suggested by Gates [7] as "one of the more useful in 
modelling the large-scale disturbances in a troposphere of constant 
e@apility 

Since it is evident that the value of 0 may approach zero or be nega- 
tive, the magnitude of @ was controlled so that the second term of (6) 


remained positive and within certain bounds. This was done by requiring 


that the lapse rate not exceed 7/8 of the dry adiabatic lapse rate, namely, 
| ee 
{o — YR (op or 
: (973 
G = (P,p a (9A) 
pp 2 
Comparison of the true O with PD Z for a few synoptic cases 
has shown that this limiting criterion is very rarely, if ever, required 


when the input data are vertically consistent. 


The second term of is approximated by the finite difference analog [8] 


aT) Aple =| SP AT + SEPA T, 


Pp Ap MP AP AYP vl (10) 


where A, and Bday are defined in figure l. 

Sigma is expressed in units of degrees Kelvin per mb. [In order to 
calculate sigma at 300 mb, temperature values at 100 mb are required. 
These data are not presently available from the Fleet Numerical Weather 
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Facility on a routine basis, and hence temperature values at 100 mb were 
approximated by either of the following schemes: 

(a) The temperature at 100 mb was set equal to the 300-mb temperature. 
This approximates mid-latitude conditions and implies that a "tropopause" 
ex1lsts Somewhere between the 300-mb and 100-mb levels, the specific press- 
ure depending on the fitted parabola, 

(b) A mean 100 mb temperature field is used based on a climatic 
summary. Since the synoptic case chosen for study was in January 1962, 
5-year mean 100-mb temperatures for January (1949-1953) were used. 

It is realized that neither of the above schemes is entirely satis- 
factory; however, the computer program iS so arranged that synoptic analy- 
ses of 100-mb temperatures can be used when available on an operational 
basis. 

In order to ut1rlize established numerical sub-routines available at 
the FNWF as well as to minimize computational time, fixed-point fractional 
arithmetic was used excluSively in the computer program. Moreover, all 
parameters had to be scaled to have magnitudes less than one. The follow- 


ing standard FNWF scaling conventions were used: 
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scaled difference equation that was programmed for the digital computer: 


A gp aciaerare aA a 
V0,-Bye2 3, a eras V Fa, >| «)- J (D, V Ix) 


where Se B..2 = “le ie A ; Px sue , ° 


(> 
Ill 


hetting i. First term on right in (12) 


( A — 
——=— 


Second term on right in (12) 


equation (12) becomes 


a | 
V x - By & =C;-Cé 


=. 5K = 1 Zeer) 
These equations were solved by an accelerated Liebmann relaxation 


iterative scheme where the (W+/ ) 1terate at the point Cam ) is given 
by 
‘Urt| S555 R~ 
is —= 
wy = aoe . (14) 
"4s 4+Bx 
Here the residual R at each Step is: 


7 ae AU 
R age VW Loe — Bx — Cr + Cyr j 


(15) 
The four pressure levels were successively scanned beginning at 


850 mb and ending with the 300-mb level, this constituting one SCAN count 
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of 4(1977 )= 7908 points in the 3-dimensional mesh. 

Considering (13) as a system of Helmholtz equations, certain condi- 
tions regarding the coefficients must be fulfilled for convergence 
('Arnason [9]) as follows. 

(a) Bx > O (16) 

mele sum of the coefficients of (x; and 774) in CAs in (ls) 


must be less than or equal to one, including round-off errors 


: A, D An | ae ; 
Bi Axe a Ae 2 (17) 


The first requirement (16) is easily met since Wk eae (see (9A)). The 
second requirement is sensitive because of the great range of the co- 
efficient Bi ,» Which can vary from a very small positive value (arbitrar- 
ily set ee to values of the order of 2000. Both of these require- 
ments are incorporated into the program. 


The convergence criterion € is defined as 


for all points of the mesh. The € is chosen so that no meteorologically 
Significant changes will occur in the omegas if iteration were continued. 
The smallest value of interest was taken as omeens ane which is approxi- 
mately equal to 0.1 cm aoe at 850 mb. 


A number of different values of € , ranging from oct to (1/6400) 


oR mb se » were tried to determine an optimum 2 he ter aad opera- 
~ 4 = | 
tional value adopted was rye = — -]0 mb sec | 


With this criterion, Cos » convergence is reached, with the synoptic 
cases studied, in 14 to 15 SCANS as defined previously. For comparison 


purposes, convergence for Cioo= a. c mb Boome is 
oO 





‘ [oO 


obtained in about 26 SCANS, 








The over-relaxation coefficient d used was determined experi- 
mentally by finding that A for which convergence was reached in the 
least number of SCANS for a given © . This optimum is very 
sharply defined and was found to be 1.414. A A differing from this 
value by roughly 10% increased the SCAN count by about 25%, and a change 


in value of only 0.1% changed the SCAN count by one. 


3. The Lower Boundary Condition 

The lower boundary condition was of special interest in devising 
this model [10]. In most previous studies, the given or calculated lower 
boundary values of omega were arbitrarily assigned to the lowest constant- 
pressure level in the model, usually 1000 mb, and then the family of 
equations was integrated with no regard to the presence or absence of 
mountainous features. The intersection of any pressure surface with the 
terrain constitutes an internal boundary within the mesh and is treated 
as such in this model. It is noted that, in order to find the inter- 
section of the terrain and a given pressure level, the pressure at the 
terrain Pr must be known. 

The input data at the various pressure levels (K) consist of Dy 


values and temperature Ty in degrees C, The definition of D is: 


= 7a fie e 


The D-+- value at the height of the terrain Zy is interpolated from the 


input data Dy by the principles of altimetry [1lj. 
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A.D =D, 7 D,. 
A, Zp = Zp.-ZLp. (18) 


MycZp = Zp. = eee 


Ds — p= Or D- = 1s = eee ° 


A.D 


Vv 


\\j 


To find A. y , the specific temperature anomaly Syy =MD/Ay Zp 
Pemeatculated. This 5;;, value is the mean for the column Pv to Pi 
Applying it to the portion from Pe to ier 

Teal - 
Substituting in (18) gives 


p. = Du=Suri( Zp. - Zr) . (19) 


, i a ee 


Having interpolated D—- , the pressure altitude is then 
Zp; = ane = 
and f— can be obtained from Zp; Since there iS a one to-one corre- 
spondence of values. 
For the last step, an interpolating polynomial is developed that 
gives the terrain pressure with accuracy within one millibar. 
The polynomial is 
ae i 4 
ia, A =p *’). (Cc: ) 20 
P y(at by") +U(c+du*) , (20) 


iS 








where a = /24,.64 C= 597. 37 
bo = 3.5562 d = 3.1272 
att 
u _ Zp, -2.7° 10° 


[See 


The scaled polynomial is then 


Pr =2° ly (44697) + 0(@4+d07)] , 


where a = 4.4229 ,1075 
Pee) 34729 alge 
€ = 3.6461.10 7 
d = 3.0539.10 
A & 
Lives A * 
A 2 
=y-Z 
y ye . 
“ Bee Pr = O25 
O,5 7220 
ya Oe 


The pressure at the terrain is calculated on a synoptic basis and is used 
in devising a marker field which controls the iterative scheme so that 
all points of the mesh which "lie" within the terrain are not relaxed 
and constitute an internal boundary. At the surface of the terrain a 
value of W,., is computed as described later. This W,,,value is assigned 
to the same grid point at all greater pressure levels, i.e., those within 
the terrain. This substitution is illustrated schematically in figure 3. 
Note that under this scheme, the Laplacian of W at (c) involves the 
uses of (Ore dir] , which generally more closely approximates the 


value of Wy than 1f zero were used. 
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4, Omega at the Lower Boundary 

The pressure at the terrain Pr is also used to calculate the omega 
at the lower boundary, W,.. The combined effect of the frictionally- 
induced vertical velocity and the vertical velocity due to the slope of 
terrain will be calculated for the volume of air bounded by the surface 
of the earth at pressure ae and at the top by a pressure eet a the 
surface frictional stress is assumed to vanish. The net inflow or con- 
vergence into this volume results in a vertical velocity Wy,through the 


upper boundary. If the equation of continuity 


- a AWWA (21) 
Op 








is integrated through the friction layer, Pt ~ Pro utd LiZinoe tine 
wind 7 
= = sere 
| 32 — | se. Stil 
V="FLoy +3 Dp Plo 7 2 Op J’ 


ee and Ty are components of the stress, there results 
































= 86 (rr-r) + $[3te- St], ws 


where i is surface stress at Pie 

and ie a ty = Cail Pie 

The first term on right in (23) is usually small in comparison with the 
other terms and is neglected. Setting Wy, -Wy SW , the 


frictionally induced vertical velocity (23) may be written as 


CO fe 


_ $) ates _ ce | . 
7 Oy Ox 


(24) 


If it is assumed (Lettau [12]) that the surface stress is opposite 
the direction of the wind, \V- , making an angle & with the isobars, 


(see figure 4), then 
° as AL Oar he AL scot AL o © 
ie a= Os COS 5 loy — lo Sing o 


The surface stress is further assumed to be proportional to the square 


of the wind at the top of friction layer in the manner of Lettau [11] 


Re 
Bodacressman [10], i.e., as, = ee Pie vz 
and Tex ~ Cp Ric U.. Vie 5 Tey =Cp Lo VioVic ; (25) 
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Then from (24) 


-— ho 
—_ B8so] dy OX zl 


where the local variation of density has been neglected. Expanding this 


8G.6) J (Cp U Mec - }O (Co vee Vie)| 


equation 


- Ny ) » -)| 
we= S84, va (iteo _ vs) ze 4, ewVix ~hh Gaye »( 26) 
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The last two terms on the right in (26) were assumed small in comparison 
with the first and were neglected. The surface friction omega can then 


beewritten as 


We = ae CoVie Foe ? fo) 


where the density has been replaced through the equation of state. 
Next, the vertical velocities due to varying terrain height will be 
considered. The value of omega due to the terrain Slope@--may be written 


iieeine form 
wr =(S2)= (SE) + VerVPr 28 


For a large 3-hourly tendency 


(<P P)_ = 5 mb per 3 br 


and will be neglected. Then 
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and omega at the lower boundary is 





Wiort Wr soe Oe, = =V,, “VP, -F EE CoVic re ° (30) 


The geostrophic wind at the lower boundary Vien is required in (30) in both 


terms. The geostrophic vorticity hie is also required and is based on Vi 6 
V.o= kx ZV D)_ (31) 
Lo £ i 


Here it is noted that Vp D must be determined at terrain pressure. Thus 
it is necessary to ''reduce" or extrapolate the De values at points sur- 
rounding the central point. These extrapolated values will be called 
(Dr), . The extrapolation will be done on altimetry principles using 
the previously calculated Sy, values for the layer in which the terrain 
lies. For example, in figure 5, at point ( érL ), the extrapolated value 


is 


(D+)p. it | — (D-) iri * (Sux ctl 





Zoi-Zprr, ir | (32) 
It follows that the gradient of D in the x-direction may be written 


(De)p int 7 (Dr) Pe a ID, Cala ap i-1| 7 Z o,d| (Sudie-(Sudéel 


(32) 
at (SuinZpr,é- a (Soles cee Need ; 


Similar considerations hold for the y-component. Both equations (32) and 


(33) are used in making calculations of Vig and eo use in equation 





(30). 
The T+ in equation (30) is also calculated by altimetry principles 


[9]. By definition, the specific temperature anomaly is | 
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S = Tate (34) 
‘Tp 


where eae the pressure temperature, is the standard-atmosphere tempera- 
ture at a given pressure altitude and T is the temperature at a given 


point. Solving (34) for T and applying to this model gives 
The pressure temperature ee given by 


Tey = Tpe- Zp, (36) 


i 
© 


where Tr 
Po 


15C is the temperature at an 


- } a 5 } 
and a = ©.5 10 °C per cm. 


Combining (35) and (36) yields 
Tr = (Tp. -a+F pz)({-Ser) : (37) 


Substituting (31), (32), (36) and the geostrophic vorticity into (30) gives 


the computational forms for W+ and O) |= 
cae. 4” 
ec = 4 (Dr)o5 Pr] > oa 


a —*Cp:|Kx V, D+ | Vv x). 
FER Clea q' Zpr)(1t+ Su.) 





(2) 


Using the finite difference analogs (7) and the standard scaling conven- 


tions (ll) gives 


Sy = Fe ol(b)p Pr] “0 


ae 





r ‘az A Dank 
2g sa? A. Be -Ang(BDp Tne 





We=- = - re (41) 
a £7 sind d?-R (Tp. -& Zp,)(ir Soe) 
The WO, ,is then LS, 6 = U3 + te - 


The computational time to obtain oe is about 40 seconds. 

The coefficient of drag Cp in (30) is the same as that defined by 
Cressman [13], and values of Cp were obtained from the Joint Numerical 
Prediction Unit, United States Weather Bureau, Suitland, Maryland. 

The height of the terrain Z7 is an area-smoothed value, and for the 
Synoptic cases presented the heights used were obtained from the FNWF. 
The terrain heights have certain deficiencies due to rather heavy smooth- 
ing. They are being replaced with values averaged over a l-degree 


latitude/longitude area, as determined by Berkofsky and Bertoni [14]. 


De Some Features of the Computer Program 

The Control Data Corporation 1604 digital computer was used in this 
investigation. It has a core memory of 32,768 words of 48 bits each, 
with a typical add time of 7.2 microseconds. 

The initial "guess" field for omega at the various pressure levels 
is obtained from the W,,, value previously computed. It was expected that 
the influence of a particular (Uat the lower boundary would decrease 
nearly exponentically with height above the terrain, hence the initial 

(c) aif 
guess field was computed by (OC eee eee ) 
where t = 1 for the first pressure level above the terrain, t = 2 for the 
second, etc. This procedure accelerates the convergence of the relaxa- 
tion scheme appreciably, as compared to an initial guess field of zero. 


Provision is made in the computer program to determine the influence 
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of the lower boundary on the omega fields in various ways. 

(a) The terrain-induced vertical velocity W-- and the frictionally- 
induced velocityWemay be calculated separately or jointly. 

(b) &,. may be set equal to zero. Here the terrain is still present 
and the lower boundary follows the variable height of the surface of the 
earth. 

(c) The computation scheme can be altered to test the effect of 
inserting the @,;,,values at an artificial lower boundary, calculated on 
the assumption that the height of the terrain is zero. The pressure at 
this lower boundary is not constant and corresponds to the sea-level 
pressure map. This simulates closely the approximations normally used 
[3], [10], [13], wherein the lower boundary values are arbitrarily 
assigned to the lowest constant-pressure level of the model. Thus, they 
are not only applied at an incorrect altitude, but ignore the presence of 
internal boundaries due to the intersection of pressure surfaces and 
terrain. 

(d) The forcing function Cz of equation (13), which involves the D- 
and T-fields responsible for the large-scale vertical velocities, can be 
set equal to zero and a solution obtained with only the lower boundary 
values Wy, present. This essentially shows the decay of the influence of 
the @),-, with height, subject here only to the static stability variations 
in the atmosphere (see equation (6)). 

Examples of the various possibilities (a) to (d) will be given in 
the next section. 

Many modes of operation of the computer program are possible and 
allow a wide variety of side products to be obtained if desired. The 
more important of these are listed below. 
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Dy at 5 levels 


qT at 6 levels 


CG at 4 levels 


Wio 

f at 4 levels 

WwW at 4 levels 
All of these can be obtained in computer-analyzed form on machine-drawn 
and labelled charts. All the charts of this investigation were produced 
in this manner except cross sections. 

It is to be noted that on the charts of Wo wherever the terrain 
projects through a pressure level, the value +lll is printed at that grid 
point, and delineates the internal boundary. 

The computer program has "smoothers" incorporated to control irregu 


larities in higher-order finite differences of the input data. These 


Smoothers are of the form 
ee 2 
Mee eet KV ALS 3 


Ma 
where A is the smoothed data, 
A is the unsmooth data; 
and KK is the smoothing coefficient and is equal to 1/8. 


Pak Cae Aas = 
This smoother is applied once to the Y De 5 V IK 4 and V ¢ (a; r) 


fields. 
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Or Discussion of Results 

The static stability 1s a bulk value, since it is derived by fitting 
a parabola to temperatures from three pressure levels of hydrostatically- 
checked analyses. The computed Gx vary as expected with regard to the 
relative stability characteristics of various air masses at lower levels 
and with respect to the proximity of the stratosphere at upper levels. 
The observed rapid variation of the static stability may not justify the 
neglect of the derivatives of Gx in certain areas. These effects are to 
be included in future computations. Some typical charts are displayed in 
Appendix C, pages 136, 139, 142, and 145. 

The influence of the lower boundary on the vertical velocity field 
was investigated, and some typical charts of the results are displayed in 
figures 6 to 10. Note that the decay of velocity with distance from the 
surface of the earth is very rapid, and even large values of w,. (at 
levels below 850 mb) do not exert a significant influence beyond the 
500-mb level. 

In many models, the vertical velocity at the lower boundary 1s 
artificially assigned to the lowest pressure level of the model. The 
present model was constrained to operate in this mode for a few tests, 
one of which is presented in figures 11 to 15. Inspection shows that, as 
expected, there is little disagreement in areas where the terrain is near 
zero height, but large discrepancies arise at those points in the vicinity 
of prominent terrain features. Areas where this disagreement is the 
largest are obviously the Cordilleras of North America and the Himalayan 
Plateau of Asia, The transparent overlay at the end of this report may 


be used on figures 6 through 10. 
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ts Vertical-Velocity Computations for a Particular Synoptic Situation 

The California snow storm of January, 1962 was selected to provide 
data for vertical-motion analyses. It was considered that this storm 
might provide results which would show the development of a situation 
which was inadequately forecast. An extensive synoptic investigation of 
this storm has been made by Ikegami [15]. 

The required input data of D values and temperatures at 1000, 850, 
750, 500 and 300 mb were obtained from FNWF, Times chosen for analyses 
were 1200Z 17 January, 1200Z 20 January, O000Z 21 January, 1200Z 21 Janu- 
ary, 0000Z 22 January, 1200Z 22 January, 1200Z 23 January, and 1200Z 
26 January. 

Appendix A contains extracted portions of hemispheric prints which 
measure 30 x 26 inches. The extracts are the same scale as the originals. 
Since various cross sections are provided, a line delineating the cross 
section path is indicated. Also, grid points three and four are under- 
lined. The Monterey Bay area lies between these two points, as shown by 
the continental outlines drawn on the chart on page 44 of Appendix A. 

This extract gives the pressure at terrain height for QOQOOZ 22 January. 
If it is desired to use the transparent overlay on these charts, shift 
the overlay to the left one grid point for proper display. 

As indicated earlier, wherever terrain projected through any pressure 
level, the value +111 was printed at that grid point. For the maps pro- 
vided, only the 850-mb level is affected. Since a point within the 
terrain was not relaxed in computing the omega fields, the terrain value 
of omega was substituted at that grid location. 

The printed values of omega multiplied by 1077 give units of mb/sec. 
Vorticity values multiplied by ome give units of eee. 
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a. Analysis for 1/ January 

(See Appendix A, pages 45 to 54 ; Appendix B, page 106) 

In general, California was under the influence of subsidence assocl 
ated with a high located to the west. However, in the northwestern United 
States a wave developed on a Stationary front oriented from northwest to 
southeast. This wave continued to develop and the associated cold front 
brought heavy rains into California. 

Rain began in the Monterey area about 1500Z 19 January and continued 
through 1300Z 20 January. Frontal passage occurred about 1000Z 20 Janu- 
ary. The pressure rise was small and no appreciable increase in pressure 
was noted for 48 hours, due to a surface trough which remained in the 
area. 

Program outputs of the region where this development occurred depict 
ascending air at all levels. Further, relative vorticity analyses show a 
maximum located northwest of the centers of rising air. 

Appendix C, page 119, shows the presence of a blocking high at the 
500-mb level. The block's orientation and the downstream cold low became 
the controlling factors in the movement of the surface wave. 

b. Analysis for 20 January 

(See Appendix A, pages 54to 59; Appendix B, page 107.) 

By 1200Z the front had passed Monterey, but intermittent showers 
continued to about 1800Z. The center of ascending air was located north- 
east of the Monterey area with higher vertical velocities at the upper 
levels. Strong subsidence to the northwest was apparently moving east- 
ward and would soon dominate the local weather. 

But, as seen in Appendix C, page 121, the previously-mentioned block 
was now fully developed. The downstream low became quasi-stationary and 
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the area of subsidence remained to the west. 
c. Analysis for 21 January 

(See Appendix A, pages 59 to 77; Appendix B, pages 108 to 110.) 

The 500-mb chart for 0000Z (Appendix C, page 128) shows essentially 
the same situation as existed 12 hours earlier. The temperature of the 
cold low was near -40C. Moreover, this tongue of cold air extended down 
to the 1000-mb level. However, the ascending air was now located south- 
east of Monterey, and higher vertical velocities were found at lower 
levels. 

The 1200Z analyses show that the center of maximum relative vorticity 
and ascending air were coincident at low levels, i.e., through 850 mb. 
Using the terrain boundary on the 850-mb chart, it appears that topography 
restrained the lower portion of the rising column, shearing it from the 
main cell. 

Rain commenced locally at 0600Z 21 January. Snow pellets began 
falling about 1400Z and snow showers began at 1615Z. Throughout the 
period ending 1000Z 22 January, snow showers, snow pellets and rain were 
frequent. 

d. Analysis for 22 January 

(See Appendix A, pages 77 to 95; Appendix B, pages 110 to 115.) 

The significant aspect at OO00Z was the change in direction of move- 
ment of the cell of rising air to the southwest. Vertical velocities 
decreased near points three and four at the 850-mb surface. Appendix C, 
pages 136, 139, and 142, shows the extent of low-level instability. How- 
ever, to the west, subsidence was clearly increasing at all levels and 
now appeared to be moving to the east. 

The 1200Z results show further movement and intensification of 
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subsidence as it proceeded into northern California. The ascending air 
continued moving southwestward. The surface low, after having reached 

a position in southeastern Nevada, moved southwest and weakened west of 
Lower California. At 1200Z 22 January, it was situated over the Los 
Angeles area. 

Qo Analysis for 23 January 

(See Appendix A, pages 95 tol00; Appendix B, page 115) 

These 1200Z charts, 24 hours later than the previous analyses, show 
points three and four were under the direct influence of descending 
motion, save for a small area of rising air due to terrain effects near 
point four. The velocity of the ascending air to the southeast generally 
decreased and the cell appeared to be weakening. 

In Appendix C, pages 146 to 156, the shift of cold temperatures and 
the trough towards the southwest and the weakening of the block at 500 mb 
are depicted. 
fer Analysis for 26 January 

(See Appendix A, pages 100 tol05; Appendix B, page 116.) 

The entire West Coast was under the influence of subsidence and 
little evidence remained of the column of rising air. Appendix C, 
page 157, shows the trough had resumed normal west-to-east movement. 

g. Conclusions 

The preceding analyses suggest that had this vertical velocity data 
been available, the forecast of fair and sunny conditions for 21 January 
would not have been issued. The vertical motion fields, together with 
Cemperature analyses, indicate the impending storm. 

The omega fields were compared to the surface analyses as trans- 
mitted by facsimile, and indicate that the direction of vertical motion 
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might have been inferred from conventional map analyses alone. However, 
the magnitude of the velocities of ascending air in the surface trough 
could hardly have been appreciated from such maps. 

Lack of time precluded computing vertical velocities for other 
Synoptic situations. Also, as interesting and encouraging as these 
indications are, it should further be pointed out that the results dis- 
played are for an area of abundant data. It remains to be seen how uSe- 
ful and accurate vertical-velocity computations are in areas of sparse 
upper-air data. 

The outputs of this program can form the basis for other investiga- 
tions in the field of numerical forecasting. Presently, the vertical 
velocities obtained here are being used to compute the transformation of 
potential energy to kinetic energy on a hemispheric scale. 

With prognostic D and temperature fields, prognostic vertical 
motions can be computed. The addition of hemispheric dew-point analyses 
and prediction would permit forecasts of cloudiness and precipitation. 
The contribution of the release of latent heat to the stability parameter 
could then be included. The obvious application of the omega fields to a 
multi-level forecast model is another possibility. 

On the hemispheric charts, numerical shading between plus and minus 
areas of vertical motion imparted an appearance of cloud patterns. 

Since facsimiles of TIROS IV cloud pictures were available, two samples 
were chosen for comparison. 

Appendix D contains the TIROS pictures and the applicable omega 
fields. Page 159is the 0501Z 17 May 1962 TIROS observation; an apparent 
vortex is indicated near 35N 163E. Here the vertical motion charts lag 


the observation by seven hours. 
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Page 165 is a 2220Z 29 May 1962 picture. A vortex is located near 
42N 100W. The omega charts are for a period about two hours later. 

The observations and representations are well correlated. Perhaps 
further study would enable the cloud type to be correlated with the 


vertical-velocity field with respect to sign, magnitude, and level. 
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APPENDIX A 


EXTRACTS FROM HEMISPHERIC PRINTS OF 


OMEGA AND RELATIVE VORTICITY FIELDS 
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CROSS SECTIONS OF OMEGA FIELDS 
PLUS TEMPERATURE, RELATIVE VORTICITY AND STABILITY 


FOR O000Z 22 JANUARY 1962 
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HEMISPHERIC CHARTS OF D VALUES 


AND TEMPERATURES PLUS STABILITY 
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